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ABSTRACT 

During 1997, the BL Lac object Mkn 501 went into an extraordinary state of high 
X-ray and TeV gamma-ray activity, lasting more than 6 months. In this paper we 
report on the TeV emission characteristics of the source in the subsequent years of 1998 
and 1999 as measured with the Stereoscopic Cherenkov Telescope System of HEGRA 
(La Palma, Canary Islands). Our observations reveal a 1998-1999 mean emission level 
at 1 TeV of 1/3 of the flux of the Crab Nebula, a factor of 10 lower than during the 
year of 1997. A dataset of 122 observations hours with the HEGRA telescope system 
makes it possible to assess for the first time the Mkn 501 TeV energy spectrum for a 
mean flux level substantially below that of the Crab Nebula with reasonable statistical 
accuracy. Excluding the data of a strong flare, we find evidence that the 1998-1999 
low-flux spectrum is substantially softer (by 0.44 ± 0.1 s t a t in. spectral index) than the 
1997 time averaged spectrum. The 500 GeV to ~10 TeV energy spectrum can well be 
described by a power law model with exponential cutoff: dN/dE oc E~ a exp (—E/Eq) 
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with a = 2.31 ± 0.22 stat , and E = 5.1(-2.3 + 7.8) sta t TeV. Within statistical accuracy, 
also a pure power law model gives an acceptable fit to the data: dN/dE oc E~ r with 
T = 2.76 ± 0.08 s tat- After presenting the 1998-1999 TeV characteristics of the source 
we discuss the implications of the results. 

Subject headings: galaxies: BL Lacertae objects: individual (Mkn 501) — galaxies: jets 
- gamma rays: observations 

1. Introduction 

During the years of its discovery as a TeV source the BL Lac object Mkn 501 showed modest 
integral fluxes, during 1995 about 1/12 of the Crab flux above 300 GeV (Quinn et al. 1996) and 
during 1996 about 1/3 of the Crab flux above 1.5 TeV (Bradbury et al. 1997). In the year 1997 
the source went into a state of surprisingly high activity in the X-ray (Pian et al. 1998) and TeV 
energy bands (Samuelson et al. 1998; Aharonian et al. 1999a; Djannati-Atai et al. 1999). During 
more than 6 months the source showed a succession of very strong flares with an average differential 
flux at 1 TeV approximately 3 times higher than the flux of the Crab Nebula. 

The extraordinary high TeV emission level as well as the significant improvement of detectors 
and analysis methods made it possible to study the temporal and spectral TeV characteristics of 
a BL Lac object with unprecedented detail. The measurements with the Stereoscopic Cherenkov 
Telescope System of HEGRA revealed that while the TeV flux varied by factors of up to 30, the 
TeV spectrum remained surprisingly stable. With a typical statistical accuracy (1 a) of the diurnal 
spectral indices between 0.1 and 0.3 no spectral variability could be ascertained with a statistical 
significance exceeding 3cr (Aharonian et al. 1999a). Furthermore, dividing the data into groups 
according to the flux level at 2 TeV did not reveal any evidence for a flux hardness correlation 
(statistical accuracy of the 1 TeV to 10 TeV spectral indices < 0.05). 

The excellent gamma-ray statistics combined with the 20% energy resolution of the HEGRA 
instrument resulted in the first detection of gamma-rays from an extragalactic source well beyond 
10 TeV, and the first high accuracy measurement of an exponential cutoff in the energy region above 
5 TeV, well into the exponential regime (Aharonian et al. 1999b). From 500 GeV to ~20 TeV the 
differential photon spectrum could be approximated by a power-law with an exponential cutoff: 

dN/dE = N (E/l TeV)- a exp (-E/E ), (1) 

with iV = (108 ± 2 stat ± 2.1 sys ) x 10~ 12 cm~ 2 s" 1 TeV" 1 , a = 1.92 ± 0.03 stat ± 0.20 sys , and E = 
6.2±0.4 stat (-1.5 + 2.9) sys TeV. 

Note that the stability of the 1997 TeV spectrum has been a matter of debate. The CAT group 
reported evidence for a flux-hardness correlation in their 1997 Mkn 501 data set: the spectrum 
of the 3 days with maximum > 250 GeV flux in their data set seemed to be harder than the 
spectrum during the rest of the 1997 observations on a statistically significant level although they 
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did not report an estimate of the systematic uncertainty on this result (Djannati-Atai et al. 1999). 
The apparent contradiction between the HEGRA and CAT results could be explained by the 
following 2 facts: (i) Since the CAT flux-hardness correlation stems from the energy range below 
1 TeV while the HEGRA constraints on the correlation are strongest at energies above 1 TeV, 
the results could indicate a stronger flux-hardness correlation below 1 TeV than above 1 TeV 
(Djannati-Atai et al. 1999); (ii) since the evidence for the flux-hardness correlation is based on the 
detection of a harder spectrum for a small number of observation nights (3 nights) it could well be 
that the flux-hardness relation found for these days does not represent the behavior of the source 
during the full 1997 flaring phase. While there are no HEGRA observations for the strongest flare in 
the CAT data sample (April 16th, 1997), the HEGRA data of the day with the second highest TeV 
emission in the CAT data sample (April 13th, 1997) also indicated a harder than average spectrum: 
the 1 TeV to 5 TeV photon index of 1.87±0.14 stat deviated by 2.7 a from the 1997 mean photon 
index of 2.25 (Aharonian et al. 1999a). Therefore a harder than average spectrum for at least 2 
of the 3 CAT nights is consistent with the HEGRA observations even if the spectrum hardened as 
much above 1 TeV as it did below 1 TeV. Note that the HEGRA observations of several later 1997 
flares with comparable and higher flux levels than that of April 13th did not show such a trend of 
spectral hardening 

In this paper we present HEGRA observations of Mkn 501 during the years 1998 and 1999 
and show how the energy spectrum evolved in the years following the 1997 flaring phase. After 
describing the telescope system as well as the dataset in Sect. 2 we report on the 1998 and 1999 
Mkn 501 light curve and energy spectrum in Sect. 3. In Sect. 4 we discuss possible implications of 
the results. 

2. The HEGRA Cherenkov Telescope System and the Data Sample 

The HEGRA collaboration operates six imaging atmospheric Cherenkov telescopes located on 
the Roque de los Muchachos on the Canary island of La Palma, at 2200 m above sea level. In 
this paper we present the 1998-1999 Mkn 501 data taken with the HEGRA Cherenkov Telescopes 
System (Konopelko et al. 1999a). The system consists of five telescopes (CT2 - CT6) which are 
operated as a single detector for the stereoscopic detection of air showers induced by primary 
gamma-rays in the atmosphere. The telescope system has been taking data since 1996, initially 
with three and later four telescopes, and since fall 1998 as a complete five-telescope system. With 
an energy threshold of 500 GeV, the telescope system achieves an angular resolution of 0.1° and an 
energy resolution of 20% for individual photons, as well as an approximate energy flux sensitivity 
v~F v at 1 TeV of 10~ n ergcm~ 2 s _1 (S/N = 5 a) for 1 hour of observation time. Note that most 
of the data presented in this paper were taken with the four telescope system. The telescope 
"CT2" was undergoing significant hardware modifications. The results presented in the following 
are based on 126 hours of data acquired between 1998, February 28 and 1999, July 7 with Mkn 501 
altitudes above 60°. The search of intraday flux variability makes use of additional 27 h of data 
with Mkn 501 altitudes between 45° and 60°. The standard data quality criteria, analysis tools, 



-4 - 



and "loose" selection cuts for spectral studies described in detail in (Aharonian et al. 1999a-b) were 
applied. 

3. Experimental results 

Fig. 1 shows the diurnal averages of the Mkn 501 integral fluxes above 1 TeV as measured 
during 1998 and 1999. Compared to 1997 when diurnal averages of up to ~ 178 x 10~ 12 cm~ 2 s _1 
(corresponding to 10 Crab units, 1 Crab unit equals 17.5 x 10 _12 cm _2 s _1 ) and a mean > 1 TeV 
integral flux of 64 x 10~ 12 cm _2 s _1 (3.7 Crab units) were observed, a striking decrease in flaring 
activity can clearly be recognized. Only 1 flare with a peak flux substantially surpassing the flux 
level of the Crab Nebula was detected on June 26/27, 1998 (MJD 50991) and June 27/28 (MJD 
50992). The mean > 1 TeV integral flux of the 1998-1999 data sample without the 1998 June flare 
was (4.8 ± 0.3 sta t) x 10~ 12 cnr 2 s _1 , a factor of 13 lower than that of 1997. 

The June 1998 flare and the correlated X-ray activity have been discussed in detail by Sam- 
bruna et al. (2000). For the day of maximum emission, June 26/27, we found evidence on the 99.4% 
confidence level that the flux increased and decreased by a factor of 2 within a time interval of 1 
hour. The spectrum of the 2 days with maximum emission seemed to be softer than that of the 
1997 observations, but the difference was not statistically significant. A fit of a power law with an 
exponential cutoff yielded the parameters (see Eq. (1)): Nq = (79± 10 sta t) x 10~ 12 cm~ 2 s _1 TeV -1 , 
a = 1.92 ± 0.30 sta t, and E = 4.1(-0.9 + 1.45) sta t TeV. For the remaining days of the 1998-1999 
observations we did not find additional evidence for intraday variability. For the typical diurnal 
observation times of about 1.5 h the 2a (S/N) flux sensitivity threshold was approximately l/6th 
of the Crab flux. For 36% (63%) of the days with observations, we found an excess with more 
than 2a (la) statistical significance. For 12% of the nights, the Mkn 501 flux was clearly below the 
sensitivity threshold for diurnal observations and an event deficit rather than an event excess was 
found in the signal region. 

We determined a "low-flux" TeV energy spectrum of Mkn 501 by analyzing 122 h of 1998-1999 
data, excluding the 4 h of data taken during the 1998-flare. Here we combine data of the 2 years 
to increase the statistical accuracy of the resulting spectrum. Analyzing the 2 years independently 
yields a very similar fitted mean flux at 1 TeV for both data sets and, within the large statistical 
errors, no significant difference in spectral shape. The 1998-1999 low-flux spectrum combined with 
the results of the fits discussed in the following are shown in Fig. 2. 

A fit of a power law model with an exponential cut-off yields the parameter values (see Eq. (1)): 
N = (10.1±1.9 stat ) x KT 12 cm -2 s- 1 TeV -1 , a = 2.31 ±0.22 stat , and E = 5.1(-2.3 + 7.8) stat TeV 
with a reduced x 2 -value of 0.71 for 13 degrees of freedom (the negative flux estimates corresponding 
to the flux upper limits in Fig. 2 are included in the fit). Within statistical errors, also a pure power 
law model 

dN/dE = N (E/l TeV)- r (2) 
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yields an acceptable fit with: N = (8.4 ± 0.5 sta t) x 10" 12 cm" 2 s" 1 TeV -1 , T = 2.76 ± 0.08, and 
a slightly larger reduced x 2 - va hie of 0.92 for 14 degrees of freedom. The systematic error on the 
absolute flux is 25%, due to uncertainty in the absolute energy scale. Systematic uncertainty on the 
shape of the spectrum is only relevant for energies below 1 TeV and corresponds to an uncertainty 
of a power law spectral index of 0.05 (see the hatched area in Fig. 2). The results of the fits to the 
Mkn 501 spectra are summarized in Table 1. 

The fit results indicate that the 1998-1999 low-flux spectrum is softer than the 1997 time 
averaged spectrum. The parameterization of Eq. (1) is not suited to assess the statistical significance 
of the spectral steepening since the fit-parameters a and Eq are strongly correlated and hence, the 
statistical errors on both parameters are large. A fit to the ratio of the two spectra shows that 
indeed the spectrum softened significantly: (dN/dE)( 1998-1999 low-flux) / (dN/dE) (1997) oc E 1 
with 7 = —0.44 ± 0.10 s tat- The systematic error A7 syst on the change in spectral index from year 
to year is smaller than the error on the shape of individual spectra since several contributions to 
the systematic error affect all measured spectra in the same way; we estimate that A7 syst < 0.05, 
i.e. it is smaller than the statistical error on 7 of 0.1. This latter conclusion is substantiated by 
observations of the persistent TeV emitter Crab Nebula (Aharonian et al. 2000). The spectral 
index determined for the two Crab observation periods 1997-1998 and 1998-1999 differs by only 
0.02±0.07 sta t- Given the statistical errors on the low-flux spectrum, it is not possible to decide 
which one of the two parameters, Eq or a, actually changed. Fixing Eq to the 1997-value of 
6.2 TeV gives a 1998-1999 photon index of a = 2.36 ± 0.10 sta t, i-e. a value 0.44 softer than that of 
the 1997 energy spectrum. Fixing a to the 1997-value of 1.92 gives a 1998-1999 high energy cutoff 
of E = (2.61 - 0.38 + 0.44) stat TeV, i.e. a value reduced by 3.6 TeV compared to the 1997 data. 

The upper panel of Fig. 3 compares the spectral energy distributions as measured dur- 
ing 1997, during the 1998-flare, and during the 1998-1999 low-flux phases. The results of the 
dN/dE oc E~ a exp (— E/Eq) fits to the data are shown by the solid lines and the fitted spec- 
tral shapes for the 1997 and 1998-flare data are compared to that of the 1998-1999 low-flux data 
by the dashed and dotted line, respectively. The lower two panels of Fig. 3 show the ratio of 
the 1998-1999 low-flux spectrum, and the 1997 and 1998-flare spectrum, respectively. While the 
1998-1999 low-flux spectrum is significantly softer than the 1997 spectrum, its shape does not 
differ significantly from the 1998-flare spectrum: the fit of a power law to the flux ratio gives: 
(dN/dE) (1998-1999 low-flux)/ (diV/d£)(1998 June flare) oc E 1 with 7 = -0.21 ± 0.12 stat . 

4. Discussion 

In this paper we present the TeV characteristics of Mkn 501 in the years of 1998-1999 following 
the major 1997 outburst phase. We assess for the first time the TeV energy spectrum of the source 
for a mean flux at 1 TeV well below the flux of the Crab Nebula. This spectrum complements the 
high-flux spectra measured so far in the sense that the Mkn 501 spectrum has now been determined 
for flux levels between 1/3 Crab units (this paper) and ~10 Crab units (Aharonian et al. 1999a). 
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Although the 1997 spectra were determined for a range of absolute fluxes differing by more than a 
factor of 5 (yet at 1 TeV by a factor of at least 2 higher compared with the flux of the 1998-1999 
low-flux data sample), we did not detect any evidence for spectral variability. In contrast we find 
evidence that the 1998-1999 low-flux spectrum is substantially softer (0.44 in spectral index) than 
the 1997 spectra. 

The Synchrotron Self-Compton (SSC) mechanism is widely believed to be responsible for the 
non-thermal X-ray and TeV gamma-ray emission in BL Lac objects (Ulrich et al. 1997). A high 
energy population of electrons embedded in a relativistic jet approaching the observer with a speed 
close to the speed of light emits X-rays as Synchrotron radiation and TeV gamma-rays as Inverse 
Compton radiation resulting from interactions of electrons with lower energy synchrotron photons. 
The recent observations of X-ray and TeV gamma-ray flares correlated to within less than half 
a day for Mkn 501 (Krawczynski et al. 2000) as well as for the other well studied TeV emitting 
BL Lac object, Mkn 421 (Maraschi et al. 1999; Takahashi et al. 1999), are naturally explained in 
the SSC scenario and strongly support this model. More detailed theoretical work is needed to 
decide whether the multiwavelength data could reasonably be described with alternative models 
assuming that the non-thermal emission of Mkn 501 is produced by hadronic interactions of a 
highly relativistic outflow which sweeps up ambient matter (Pohl &: Schlickeiser 2000), by inter- 
actions of high energy protons with gas clouds moving across the jet (Dar & Laor 1997), or by 
interactions of extremely high energy protons with ambient photons (Mannheim 1998), with the 
magnetic field (Aharonian 2000), or with both (Miicke & Protheroe 2000). Note that in the proton 
synchrotron model the stable spectral shape is explained by the self-regulated synchrotron cutoff, 
while steepening in the low-flux state could be explained by the drop of the acceleration rate of 
protons (Aharonian 2000). While the interpretation of TeV gamma-ray spectra is hampered by the 
unknown modification of the spectrum due to intergalactic extinction as discussed, e.g. by Aharo- 
nian et al. (1999b), the temporal evolution of the gamma-ray flux and the spectral shape is free of 
this uncertainty and should be explained by models of the origin of the TeV radiation. 

Assuming for the time being that the mechanism responsible for the X-ray and TeV gamma-ray 
emission has been identified, a next step concerns the understanding of how the emission region(s) 
is (are) embedded in the jet, how an emission region evolves with time, and where the energy 
which is ultimately converted into the observed non-thermal X-ray and gamma-ray radiation comes 
from. Within SSC models the stability of the TeV energy spectrum during the 1997 flares can be 
explained by a spectrum of accelerated electrons (and possibly positrons) which is stable throughout 
the whole flaring phase (Aharonian et al. 1999a; Konopelko et al. 1999b; Krawczynski et al. 2000). 
The steepening of the TeV energy spectrum reported in this paper can be accounted for, by e.g., 
(i) a steepening of the spectrum of accelerated particles streaming into the emission region; (ii) the 
shift of the break in the electron spectrum (caused by the synchrotron and Inverse Compton cooling 
of the electrons, or by a lower energy cutoff of the spectrum of accelerated particles) towards lower 
energies; or, by (iii) a shift of the maximum energy of accelerated particles towards lower energies. 
Detailed modeling of the now available very detailed multiwavelength data on Mkn 501 as given 
e.g. by Pian et al. (1998), Djannati-Atai et al. (1999), Krawczynski et al. (2000), and Sambruna et 
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al. (2000) should make it possible to identify the origin of the spectral curvature observed at X-ray 
and TeV energies during 1997 and to determine which properties of the emission region(s) changed 
from 1997 to 1998-1999. 
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Fig. 1. — The 1998 and 1999 light curve of Mkn 501 as measured with the HEGRA telescope system. The 
dashed line indicates the steady emission level of the Crab Nebula. Each data point is a diurnal average. 
Note that during 1997 the source reached > 1 TeV integral flux levels of about 10 Crab units: HEGRA 
detected a maximum flux of (178 ± 10 sta t) x 10~ 12 cm" 2 s" 1 on June 26th, 1997. The 1997 mean > 1 TeV 
flux level was 64 x 10~ 12 cm~ 2 s _1 (3.7 Crab units). Upper limits are given at the 2 a confidence level. 
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Fig. 2. — The figure shows the 1998-1999 low-flux spectrum measured with the HEGRA telescope system. 
The systematic error on the curvature of the spectrum is given by the hatched region. The additional 
uncertainty of 15% in the absolute energy scale is not shown here. The solid and dashed lines give the fit 
results of a pure power law and a power law plus exponential cutoff model, respectively (see text). Upper 
limits are given at the 2 o confidence level. 
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Fig. 3. — The upper panel shows the spectral energy distributions measured during 1997, the 1998-flare, 
and the 1998-1999 low-flux period. The solid lines show the fit results dN/dE oc E~ a exp (-E/E Q ) with 
the pairs {a,E ) of (1.9,6.2 TeV) for 1997, (1.9,4.0 TeV) for 1998-flare, and (2.3,5.1 TeV) for the 1998-1999 
low-flux data sample. The dashed and dotted lines show the shape of the 1997 and 1998-flare spectra overlaid 
on the 1998-1999 low-flux spectrum. The lower 2 panels show the ratios r(E) of the 1998-1999 low-flux and 
the 1997 and the 1998-flare spectrum, respectively, as well as the results of power law fits r(E) oc E 1 . All 
upper limits are given at the 2 a confidence level. 
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Table 1. Results of Mkn 501 spectral fits. Statistical errors only, see text for systematic errors. 



Data Sample 


N a ' b 


r a 




E a 


Xr/ d ofs 


1997 


108 ±2 




1.92 ±0.03 


6-2t?J 


1.66/14 


June 1998 Flare 


79 ± 10 




1.92 ±0.30 




0.54/13 


1998-1999 Low-Flux 


8.4 ±0.5 


2.76 ±0.08 






0.92/14 




10.1 ± 1.9 




2.31 ±0.22 


s.ity 


0.71/13 



a Model parameters of Eqs. (1) and (2) 

b Flux normalization constant in (10~ 12 cm~ 2 s _1 TeV -1 ) 



